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We discuss the possibility of new sources of flavour violation in the squark sector of supersymmetric 
models in the context of the dark matter relic density. We show that the corresponding non-minimal 
flavour violation terms in the squark mass matrices can have an important impact on the thermally 
, averaged (co)annihilation cross section of the neutralino, and in consequence can modify its predicted 

^ relic density. We discuss in detail the relevant effects and present a numerical study of neutralino 

I ' annihilation and coannihilation in this context. We also comment on the LHC phenomenology of 

the corresponding scenarios. 
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I. INTRODUCTION 
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' Among the numerous extensions of the standard model of particle physics, supersymmetry ranks among the most 
, popular ones. In particular, the Minimal Supersymmetric Standard Model (MSSM) is probably the best studied 
• ' scenario of new physics. It allows to cure the hierarchy problem by stabilizing the Higgs mass and leads to gauge 
coupling unification. Moreover, it includes promising candidates for dark matter, whose presence remains the most 
compelling observational evidence for physics beyond the standard model. 

Nevertheless, several open questions remain, e.g., concerning the flavour structure of the theory. While models with 
minimal flavour violation (MFV) [iHsl assume that the mechanism of flavour violation is the same as in the standard 
model, the framework of non-minimal flavour violation (NMFV) allows for new sources of flavour mixing, depending 
on the exact mechanism of supersymmetry breaking. In the former case, the rotation of the Yukawa couplings from 
; ^ ' gauge to mass eigenstates remains the only source of flavour violation, and thus all flavour- violating interactions are 
' parameterized through the CKM- and PMNS-matrices as in the standard model. For NMFV, the terms originating 
from the additional sources are not related to these matrices, such that they are considered as additional parameters 
at the SUSY scale. 

In recent years, supersymmetric scenarios beyond minimal flavour violation have received considerable attention 
in the community, especially in the context of signatures at current or future colliders. Concerning (s)quark flavour 
violation, the production and subsequent decays of squarks and gluinos at the Large Hadron Collider (LHC) have 
been studied, e.g., in Refs. ^4-9]. Apart from the production of superpartners at colliders, the flavour- violating terms 
also appear in the (co)annihilation cross section of the neutralino, which is needed in the calculation of its relic density 
for a given scenario. 

In case of neutralino pair annihilation into quarks, the squarks appear as internal propagators. Additional flavour- 
violating terms can then increase the relative contributions of these diagrams since the mass splitting of the squarks 
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is modified. Moreover, NMFV allows for efficient annihilation into final states, that are forbidden in the case of 
MFV. Flavour violating effects are also important in the case of coannihilations of a neutralino with a squark, since 
the latter is then an external particle. The importance of such processes crucially depends on the mass difference 
of neutralino and squark. The increased mass splitting of the squarks can therefore have an important impact on 
coannihilation processes. Finally, also in this case new final states are opened, leading to additional coannihilation 
channels. Recently, the impact of non-minimal flavour violation in the sector of sleptons on the coannihilation of a 
neutralino with a slepton has been discussed in Ref. fio!]. 

The aim of the present paper is to provide a study of quark flavour violation in the context of neutralino dark 
matter. In this context, possible flavour-mixing effects are generally not considered in the literature. We present a 
detailed analysis of neutralino pair annihilation and neutralino-squark coannihilation in the MSSM beyond MFV. In 
Sec. ini we will briefly introduce the MSSM with NMFV in the sector of squarks and discuss its parameterization. 
The role of generation mixing in the context of neutralino (co)annihilations is discussed in detail in Sec. IIIIl Sec. IIVI 
is then devoted to numerical examples in the context of neutralino (co)annihilation and its relic density. A discussion 
of LHC phenomenology for the corresponding scenarios follows in Sec. |Vl Finally, conclusions are given in Sec. IVII 



II. THE MSSM BEYOND MINIMAL FLAVOUR VIOLATION 



In the standard model, the only source of flavour violation are the Yukawa interactions, since their diagonalization 
leads to a mismatch between flavour and mass eigenstates of quarks and Icptons. The flavour structure of the quark 
sector is very well described by the Cabibbo-Kobayashi-Maskawa (CKM) matrix, which only appears in charged 
currents, while flavour changing neutral currents are strongly suppressed. In supersymmetric theories with minimal 
flavour violation (MFV), the Yukawa matrices remain the only source of flavour violation, so that all flavour violating 
interactions of squarks are also related to the CKM-matrix. However, new sources of flavour violation may be present 
in supersymmetric models, especially if they are embedded in a grand unification framework. Depending on the exact 
realization and the involved representations, specific relations to the Yukawa matrices can lead to fiavour non-diagonal 
entries in the soft-breaking terms. These are not related to the CKM-matrix and the corresponding framework is in 
consequence referred to as non-minimal flavour violation (NMFV). 

Considering the most general flavour structure, the squark mass matrices at the electroweak scale take the form 

M'^ ('^■'■'■'■■^'^'■'^] (1) 

for q = u,d, respectively. Their diagonal blocks are given by 

■^Irr = Ml + ml + edmlsm'^9wCOs2l3, (2) 

A^Ill M? +m^ + m| cos 2/3(/d- erf sin^ 6*11/), (3) 

■^l,RR = Af? + + e^rnl sin^ 6*^ cos 2/3, (4) 

-^l,LL = ^ckmM? Vftj^M + ™« + cos2/3(/„ - e„ sin^ Ow), (5) 

where Mq, Mfy, and Mjj are the soft-breaking mass terms of the squarks. The diagonal mass matrices of up- and 
down-type quarks are denoted m„ and m^. Due to the SU(2) symmetry, the left-left entries are related through the 
CKM-matrix Vckm- The above expressions also involve the mass mz of the Z-boson, the fractional electric charge 
Cq and the weak isospin Ig of the (s)quark, the weak mixing angle 6w, and the Higgs-mixing parameter /3 defined 
through the ratio of the vacuum expectation values of the two Higgs doublets, tan/3 = Vu/vd- 
The off-diagonal blocks of the matrix in Eq. ((TJ are given by 

■MIrl = (A^Ilr)^ = ^Tu^^i*m^ cot 13, (6) 

-^Irl - (-^Ilr)' - ^T^-M*mrftan/3, (7) 

where /i is the Higgs mass parameter. The trilinear matrices Tij,d are related to the soft-breaking matrices A^^d 
and the respective Yukawa matrices Yu.d through {Tij^d)^^ = (^u,ti)y {Yu.d)ij- At the GUT scale, the usual CMSSM 
condition (A^)^^ = (^d)33 = applies, and the numerical values for Au,d at the SUSY scale are obtained through 
renormalization group running. All parameters appearing in Eqs. ([2]) to ([7]) are understood to be in the super-CKM 



3 



basis, where the neutral currents are flavour-diagonal and the quark (but not the squark) fields are in the mass 
eigenstate basis [ll|, [l^l ■ 

In order to have a scenario-independent and dimensionless parameterization of flavour-mixing, the ofF-diagonal 
entries are usually normalized to the diagonal ones according to 



= (10) 

- f,iTD),/^{M^)A^ll),,, (12) 

- ^(4),/v/K).(^).' (13) 

- ^^Tl)J^iH)A%),r (14) 



V2^ 



The normalization factor is defined in terms of the corresponding diagonal elements of the soft-breaking matrices. We 
emphasize that the following numerical analysis is based on the diagonalisation of the full 6x6 mass matrices. This 
is realized by introducing two rotation matrices, such that 

n^Mlnl = diag(m?^,...,m?J (15) 

with the mass order rriq-^ < ■ ■ < nigg for q — u,d, respectively. The rotation matrices appear in the couplings of 
squarks with other particles, and, in consequence, the flavour-violating elements will influence observables like decay 
widths or production and annihilation cross sections. Analytical expressions for couplings including squark generation 
mixing can, e.g., be found in Refs. i4-6]. We shall discuss the relevant couplings for our analysis in more detail in 
Sec. ini 

A large variety of experimental measurements puts constraints on the parameter space of new physics models. 
Below are summarized all the constraints that will be considered (at the 95% confidence level) in this study. The 
most important one for study this is naturally the relic density of cold dark matter. Combining data from the WMAP 
satellite and other cosmological measurements, the relic density of dark matter in the universe is constrained to [isj 

f^CDM/i^ = 0.1126 ±0.0036, (16) 

where h denotes the present Hubble expansion rate Hq in units of 100 km s^^ Mpc^^. 

Then, searches for superpartners at LEP and Tevatron lead to the following mass limits for Higgs bosons, neutralinos, 
charginos, squarks, and gluinos: m^^o > 114.4 GeV, m^o > 46 GeV, m~± > 94 GeV, TOj^ > 96 GeV, rrig > 308 GeV 
[l^ . Moreover, recent results from the Large Hadron Collider (LHC) lead to more stringent limits for squarks and 
gluinos within the constrained MSSM [l^. These limits are based on the hypothesis of minimal fiavour violation, so 
that we do not take them into account explicitly in the present study. Note, however, that the scenarios considered 
in the following feature rather heavy gluinos. For the lightest Higgs boson, we require irifio > 111.4 GeV, taking into 
account a theoretical uncertainty of 3 GeV [16,] . 

Moreover, precision measurements in the sector of D-, B-, and K-mesons constrain some of the flavour-violating 
elements in the mass matrices. In particular, flavour mixing involving the first generation of squarks is severely limited 
[m, [13, ■ We therefore focus on flavour mixing between the second and third generation squarks. 

The most relevant constraints on such flavour mixing are listed in Tab. |T] together with the current experimental 
measurements and the theoretical error estimate. If not indicated otherwise, they are taken from Refs. [3, [l8| . 
They include branching ratios of rare decays, B-meson oscillation measurements, the electroweak p-parameter, and 
the anomalous magnetic moment of the muon. For the latter, taking into account recent calculations which bring 
the standard model theoretical expectation closer to the experimental measured value p^ . we use only the upper 
bound given in [l^ as a constraint. In the following study, the most important limits are imposed through the precise 
measurements of the rare decay & — ^ 57 and the B-meson oscillation parameter AMb^ ■ 
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TABLE I: Experimental constraints on the MSSM parameter space, in particular on quark flavour violating elements. 



10* X BR(6 sj) 
10* X BR(B, ^ M+M") 
AMb, [ps-i] 

Ap 
10" X Aa„ 



Exp. value Exp. error Theor. uncertainty 



3.55 


±0.26 


< 5.6 [21J 




17.77 


±0.12 


< 0.0012 [23] 




255 


±80 



±0.23 [20] 
±3.3 [22] 





FIG. 1: Feynman diagrams for the annihilation of neutralinos into fermion pairs through the exchange of a neutral Higgs boson 
Hi — h'^,H'',A'' (left), a Z^-boson (centre), or a sfermion (right). The corresponding ii-channel diagram, obtained through 
crossing, is not shown. 

III. IMPACT ON THE RELIC DENSITY OF DARK MATTER 

The relic abundance in our universe of a stable dark matter candidate can be evaluated by solving the Boltzmann 
equation 



dn 
'dt 



3Hn 



iw) (r 



(17) 



where n is the number density of the relic particle, H the (time-dependent) Hubble expansion rate, and ricq the number 
density in thermal equilibrium. All information concerning the particle physics model parameters is contained in the 
annihilation cross section tTann multiplied vi^ith the relative velocity v of the annihilating particles. This product has to 
be convolved with the velocity distribution of the non-relativistic dark matter particle in order to obtain the thermally 
averaged cross section ((Tann'f^)- 

Denoting the mass of the dark matter candidate by toq and taking into account a set of N potentially co-annihilating 
particles with masses rrii (i = 1, . . . ,N) such that toq ^ "^i < • • ■ < ttin, the thermally averaged annihilation cross 
section can be written as [2J, 



N 



cq cq 



"cq 



3/2 



exp 



(ruj + nij — 2mo) 
T 



(18) 



Here, denotes the equilibrium density of the particle i and i — refers to the dark matter candidate. The cross 
sections tiy relate to the different coannihilation processes within the ensemble of particles and Vij is the relative 
velocity between the particles i and j. Moreover, gi denotes the number of degrees of freedom of particle i and gcs is 
a normalization factor. 

From Eq. (|18p it becomes immediately clear that the mass differences between the annihilating particles play a 
crucial role. Due to the exponential suppression, the coannihilation between two given particles i and j will only lead 
to a significant contribution, if the two masses rui and rrij are nearly degenerate [241 ■ 

In the following discussion, we assume that the lightest neutralino is the lightest supersymmetric particle (LSP) 
and therefore the dark matter candidate. In wide regions of the MSSM parameter space, the pair annihilation of two 
neutralinos into standard model particles is the dominant process. The diagrams for annihilation into quarks, i.e. 
where flavour violation in the (s)quark sector can become relevant, are shown in Fig. [T] At the tree-level, squarks can 
then appear only in internal propagators in case of annihilation into quark-antiquark pairs, i.e. XiXi ~^ through 
the exchange of a squark in the t- or u-channel [IE 113] ■ 

Going beyond minimal flavour violation, the mass splitting of the involved squarks is increased due to the additional 
off-diagonal elements in the mass matrix. In particular, the lightest squark mass eigenstate (purely stop-like in the 
CMSSM with MFV) becomes lighter with increasing flavour mixing. Its contributions to neutralino pair annihilation 
through t- or u-channel exchange are therefore enhanced. 
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FIG. 2: Feynman diagrams for the coannihilation of a neutralino with a sfermion into a fermion together with a Higgs boson 
H — hP,H^\A^,H^ (top) or a gauge boson V — j, Z'^ , g (bottom). These processes proceed through the exchange of a 
fermion (left), a sfermion (centre) or a gaugino (right). The u-channel diagram is not present in case of the gluonic final state. 



Apart from the impact on the squark mass eigenvalues, the flavour- violating terms discussed in Sec. |TT] directly affect 
the neutralino-squark-quark coupling, which is present in the t- or u-channel diagram. The analytical expressions for 
the left- and right-handed parts of this coupling are given by Q 



(e„ - T„) sin OwMii + r„ cos Ov/M-a 

(ed - Td) sin 6*1^^1 + 7d cos9w^f^2 
TO„, cos 9w 



TJU* , """fc """" " '''' Kf -TJU 



^d* , rud^ cos Ow -nd* 
^^^^ 2mwsinr'' 



■ a \r -nd . "^d^ COsOw 



2mvKsin/3-^^'^^''=+3)' 



(19) 
(20) 
(21) 
(22) 



with the same notations as in Sec.|TTl Flavour mixing effects arise through the squark rotation matrix TZ"^ (q — u,d). 
This can allow for new annihilation channels, that are closed in the case of minimal flavour violation. Such channels 
can, e.g., be xx — >■ cc through exchange of a squark ui which is now a mixture of c and t. In the case of MFV, this final 
state is only possible through exchange of a heavier c and therefore suppressed. Another example is annihilation into 
a mixed final state, xX ~^ ct, which is forbidden in MFV. The discussed enhancements and new channels increase the 
total annihilation cross section, which in turn decreases the predicted relic density of the neutralino. The diagrams 
with s-channel exchange of a Higgs or gauge boson remain insensitive to squark flavour mixing. 

Let us now turn to the case of neutralino-squark coannihilation. The possible final states are a quark together 
with a Higgs or a gauge boson. The relevant Feynman diagrams at the tree-level are depicted in Fig. [5J The main 
impact from non-minimal flavour violation will be through the modified squark mass spectrum. As already stated 
above, the mass difference between neutralino and the lightest squark enters the calculation of the corresponding 
coannihilation cross section exponentially. When the squark mass approaches the neutralino mass due to increasing 
flavour mixing, this can significantly enhance the contribution from the corresponding coannihilation with respect to 
the case of minimal flavour violation. 

Again, also the flavour-violating couplings can have subdominant effects on the coannihilation processes. Each of 
the diagrams depicted in Fig. [5] contains the squark-quark-neutralino coupling already discussed above. Moreover, 
the couplings of squarks to Higgs- and massive gauge bosons are sensitive to flavour-violating effects. In the mass 
eigenstate basis, the couplings of squarks to a Z°-boson are given by [6| 



C 



92 



COS 6w 



^I'^lj'^tk - eg sin^ OwSjk 



(23) 



for q ~ u,d. Here, pj and pk denote the momentum of qj and qt, respectively. The interactions of squarks with a 
photon or a gluon are flavour-diagonal and are therefore not discussed in detail here. 
The couplings of two up-type squarks with the light scalar Higgs boson are given by @ 



h°UjUk 



92 

2mw 



E 



sin(« -f /3) (1 - i tan^ 0w)n%nt* + ^ tan^ ^?w7^"(,+3^7^; 



j(«+3)'^fc(i+3) 



6 



+2 



sin/3 



2 'T? ti* 



"fc(i+3) 



sin/3 



''"shT^ [^i(»+3) (^c^)*' T^ki + niii^^y 



(24) 



From this expression, the coupling to the heavy scalar Higgs is obtained through the replacements /i" — >■ and 
a — a + 7r/2. Moreover, couplings of down-type squarks to the neutral scalar Higgses are obtained by replacing 
Ui — ^ di and sin/3 cos/3. Finally, the couplings of up-type squarks to a pseudoscalar Higgs-boson are given by 0] 



C 



AO 



■ 92 



j(«+3)' 



3 

^.,7^f * + cot /3^ ^ 7ef(,+3) (Tc;),,7e|* + h.c 



(25) 



Again, the expressions for down-type squarks can easily be obtained through Ui — >■ di and cot /3 — >■ tan /3. 

The effects of the modified mass eigenvalues and the modified couplings are superimposed. Since the two effects are 
linked together through their same origin (see Eq. ([T5|)). the separate impacts on the (co)annihilation cross section 
and the neutralino relic density cannot be disentangled. However, some general features can be expected. The effect 
of the modified squark mass eigenvalues on coannihilation is expected to be stronger than in the case of neutralino pair 
annihilation due to the exponential factor already mentioned above. Moreover, the squark is here an external particle, 
and the impact of its mass on the phase space is more important than the mass in the t- or w-channel propagator. 

The impact of the modified fiavour contents of the involved squarks, i.e. the effect of the rotation matrix in the 
coupling, is expected to be smaller than the mass effect. This is again due to the exponential factor in Eq. (ITS)) . Note 
also that, the mixing being unitary, the newly opened channels can be (partially) compensated by the simultaneous 
diminution of other contributions. The compensating contribution can, however, turn out to be forbidden in specific 
kinematical configurations and the impact of the new contributions can be significant. This is in particular the case 
when the neutralino is too light to annihilate into top-quark pairs, i.e. for to^o < rrit. The fiavour violating elements 
lead then to a c admixture in the lightest squark, which then allows for neutralino pair annihilation into top and 
charm quarks. 

Note that there can also be coannihilation of a neutralino with an up- (down-) type squark into a charged Higgs boson 
or a W-boson together with a down- (up-) type quark. In this case, the u-channel diagram includes a chargino 
propagator and in consequence the corresponding chargino-squark-quark coupling, while the s- and t-channel diagrams 
involve couplings of up- and down-type squarks to the charged Higgs or W-boson. Analytical expressions for these 
couplings can be found in Refs. Since they are rather similar (with obvious replacements, e.g., concerning gaugino 

mixing) to the interactions given in Eqs. (|T^ to (P^ . they are not displayed in detail here. Note, however, that these 
couplings explicitly depend on the CKM-matrix. The general argumentation given above remains unchanged. 



IV. NUMERICAL ANALYSIS 



The following numerical analyses are mainly based on the constrained MSSM with the five parameters mo, toi/2, 
Aq, tan/3, and sgn(^). We also consider variants of this model featuring non- universal Higgs or gaugino masses. 
Starting from the high-scale parameters, the soft-breaking terms at the scale Q = 1 TeV [2g| are obtained through 
renormalization group running using the public program SPheno 3 [16]. At the same scale, we introduce the non- 
diagonal entries in the squark mass matrices as discussed in Sec. |lTl The physical mass spectrum is then calculated 
again using SPheno, which takes into account the general flavour structure. The same code is also used for the 
evaluation of the constraining observables mentioned in Sec. [Hi again taking into account squark generation mixing. 
For the standard model parameters, we refer the reader to Ref. [14j. The pole mass of the top-quark is taken to be 
TOtop = 173.1 GeV according to recent measurements from DO and CDF [23|. The CKM-matrix is taken in the usual 
Wolfenstein parametrization with the recent values A = 0.2253, A = 0.808, p — 0.132, and fj = 0.341 [T^ . 

Making use of the SUSY Les Houches Accord [s^l , the mass spectrum and related mixing parameters are transferred 
to the public program micrDMEGAs 2.4 [SlJ in order to evaluate the relic density of the neutralino. The calculation 
of the annihilation cross section is done by the program CalcHEP [s^], where we have implemented the MSSM with 
squark generation mixing as discussed in Sec. [Ill The corresponding model files have been obtained using the package 
SARAH [33|. We also include important effects from the running strong coupling constant and running quark masses, 
as they are also included in the default implementation of the MSSM in micrOMEGAs / CalcHEP. 
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Ao = - 500, tan p = 1 0, n > 0, 5„„r 23 = = - 500, tan p = 1 0, n > 0, 8^„„ 23 = 0.98 




Mill i^^^^^M^^P^^^^^^^^^^^^^^^^I U ^ ^ I I I U 
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FIG. 3: Cosmologically favoured region and related exclusion limits in the (mo, mi/2) plane of the CMSSM for S^^i = (left) 
and 52;^^^^ = 0.98 (right). 

A. Constrained MSSM 

In order to illustrate the numerical influence of flavour violating elements, we start by analyzing the neutralino relic 
density within the constrained MSSM (CMSSM), where we allow for flavour violation between the second and third 
generation of up- type squarks in the right-right chiral sector. In Fig. [3l we show typical scans of the mo-mi/2 plane 
for fixed values of Aq — —500 GeV and tan/3 — 10 and for positive values of /i. The cosmologically favoured region 
of parameter space according to Eq. (|16p together with the relevant constraints discussed in Sec.|TT]are shown for the 
case of minimal flavour violation (MFV, 523^^ = 0) and for the case of important off-diagonal elements, ^23^^ = 0.98. 

In the case of MFV, the most stringent constraints on this parameter plane are due to a charged dark matter 
candidate (low mo), tachyonic solutions of the renormalization group equations (high toq and low mi/2) as well as 
the constraints from b ^ S"/ and the lightest Higgs mass (low mass region). The cosmologically favoured region of 
parameter space is divided into several distinct regions: the so-called focus point region (high mo, not visible here), 
the resonance of the light Higgs boson (low mi/2 and moderate mo), and the coannihilation region (close to the 
exclusion due to a charged dark matter candidate), where the neutralino mass is close to the stau mass. 

In the corresponding figure for the NMFV-case, we depict the same constraints together with the relative contribu- 
tion from new (co) annihilation channels as discussed in Sec. IIIIl In this case, this involves neutralino pair annihilation 
into a mixed charm-top final state and coannihilation of a neutralino with the lightest squark ui. 

In the latter corresponding region (mj^/2 ^ 450 GeV), where the relic density constraint is fulfilled, the mass 
difference between the lightest squark and the neutralino is about 30 GeV, as can be seen from the left panel of Fig. 
m where we show the cosmologically favoured regions of parameter space in the plane of the physical masses. The 
dominant annihilation processes are then x^^i — >■ gt (30%) and uiUi — )■ gg (25%). Two other important processes 
arc neutralino annihilation into pairs of top quarks (10%), and Xi^i — > gc (15%). Note that the presence of a charm 
quark in the final state is a genuine effect of flavour violation. Indeed, as a consequence of the off-diagonal elements 
in squark mass matrices, the lightest up-type squark is here a mixing of tj^ and cn (with a small admixture of i^), 
opening up the (co)annihilation into charm-quarks. 

For lower masses (e.g. mo ~ 200 GeV and mi/2 ^ 400 GeV), coannihilation processes such as Xi^i ~^ 9t/c are still 
important (20%). However, the squark being much lighter (m^^j ~ 190 GeV), the squark pair annihilation uiUi — >■ gg 
is now subdominant. Moreover, the neutralino mass of m^o 160 GeV (see Fig. |4] left) forbids annihilation into 
top quark pairs. As a consequence, the flavour violating process XiXi ~^ tc{ct), which is kinematically allowed and 
enhanced by the rather light squark in the i-channel propagator, becomes important (40%). This is represented by 
the green area in the left part of the plot. Notice the cut at mi/2 ~ 420 GeV, which corresponds to m^o w mj. 
For m^o > mt, neutralino annihilation into top quark pairs is kinematically allowed, and the tc(cf) final state is 
suppressed. This can also be seen in relation to the physical neutralino and squark masses in Fig.[5]left. For low mi/2 
but large mo, the squark being heavier, coannihilation is not relevant and neutralino annihilation into tc{ci) is less 
important. Therefore, even if the relative contribution of this channel is still important, its absolute contribution is 
not large enough to satisfy the relic density constraint. 
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FIG. 4: Cosmologically favoured region and related exclusion limits for — 0.98 in the (m^o, — m^o) plane for fixed 

Ao = -500 GeV and tan^ = 10 (left) and in the {S^i^^, Ao) plane for fixed mo = 200 GeV and mi/2 = 400 GeV (right). 

In the region excluded by BR(5 — )■ sj) most of the deviation from the standard model value comes from large negative 
chargino contributions due to the smallness of the stop and/or chargino mass. There is, however, no significant effect 
coming from the flavour violating parameter , since BR(6 57) constrains mainly flavour violation in the 

left-left sector. 

Let us now discuss the interplay of helicity mixing and additional flavour mixing. The former is induced through 
the trilinear matrices Tu (see Eq. ([6])) and thus the GUT-scale parameter Aq, while the latter is included at the 
electroweak scale through the parameter ^23^^- case of MFV, i.e. for $2^^ = 0, a rather large |^o| is needed 

in order to decrease the stop mass close to the neutralino mass, and therefore allow for efficient coannihilation. For 
sizeable additional flavour mixing, the coannihilation is important already for lower values of Aq, since the squark 
mass splitting is then increased by the off-diagonal elements in the mass matrix. 

This is illustrated in the right graph of Fig. 21 where the constraints, cosmologically favoured regions, and different 
contributions to the annihilation cross section are shown in the {Ao,52^^) plane. The mass splitting of the squarks 
depends strongly on both of these parameters, which therefore have a competitive effect on the light stop mass. 
As a consequence, as explained above, one of these parameters has to be large in order to allow for an important 
coannihilation contribution. On the other hand, the flavour violating effects are only related to 623^^ ■ Therefore the 
flavour violating neutralino annihilation processes depend mainly on this parameter. The only possibility to satisfy 
simultaneously the relic density and BR(6 — >■ 57) constraints is for very large (^23^^ ^ rather low ^o- This is 
explained by the strong dependance of BR(6 — >■ 57) on the squark mass spectrum, and therefore on Aq. Contrary, and 
as explained above, BR(6 — 57) does not depend on any flavour mixing among right up-type squarks, and the mass 
effects become important only for very large values of 62^ . It has been checked that the other constraints described 
in Tab. U are fulfilled in the whole parameter space shown in Figs. |3] and 21 Moreover the calculated spectrum is 
compatible with the mass limits given in Sec. [Ill except for the stop in some regions where it is the LSP (i.e. already 
excluded) . 

Next, we study the possibility that not only the parameter 82^^ is large, i.e. of 0(1), while all others are small, 
which might not be very natural. We therefore show in Fig.Othe cosmologically favoured region and related exclusion 
limits in the {^2^'^ ,^23"^) plane for fixed toq — 200 GeV, mi/2 = 400 GeV, and Af) = —500 GeV. We observe that 
the second flavour- violating parameter 82^^ can reach values up to 0.15 before being constrained by the lower Higgs 
mass bound of 111.4 GeV. Similarly, the RL and LL parameters (not shown) are restricted by the FCNC process 
6 — >■ 57 to values below 0.15 and 0.1, respectively, as would be the LR parameter if one applied this limit at the two 
(not three) sigma level. 

In Fig. [51 we show for a given parameter point the neutralino relic density and the contributing processes as a 
function of the flavour- violation parameter Sl^^^- While for the case of MFV, this scenario is cosmologically strongly 
disfavoured with fl^oh'^ ^ 20, the relic density decreases with increasing flavour mixing to reach the favoured value 

of ri^o/i^ 0.11 for (^23^^ ^ 0.98. For low values of 62^^, the annihilation is dominated by lepton final states (about 
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FIG. 5: Cosmologically favoured region and related exclusion limits in the (<523^^, ^23^^) pla-ne for fixed mo = 200 GeV, 
mi/2 = 400 GeV, and Aq = -500 GeV. 
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FIG. 6: Relic density of the neutralino (left) and contributing processes (right) as a function of ^23 for mo = 1500 GeV, 
mi/2 = 680 GeV, Aq = -500 GeV, tan/3 = 10, and ^ > 0. 



75%), which do, however, not lead to a sufficiently enhanced annihilation cross section. The subleading channel is 
annihilation into top-quark pairs (about 25%). For S^^ > 0.2, flavour violation effects start to manifest by opening 
the channel XiXi ~^ elite). The relative contribution of this process amounts to almost 40% at S^^^ ~ 0.8. For 



^M^RR ^ g annihilation into top-quarks is significantly enhanced due to the lighter squark in the ^-channel 

propagator, so that this channel remains more important than the newly opened annihilation into top- and charm- 
quarks. All contributions from neutralino pair annihilation drop at ^23^^ Q-Sib when the squark ui becomes light 
enough for efficient coannihilation. The corresponding total relative contribution amounts to about 60%. When the 
squark becomes even lighter, also squark pair annihilation into gluon pairs plays an important role (see Eq. (jl8p ). 
leading to relative contributions of about 90% at most. 

For this discussed scenario, the favoured relic density of the neutralino is achieved through important coannihilation 
for rather large values of the flavour mixing parameter S^^^- Note that, depending on the exact parameter point 
under consideration and the corresponding relic density in the MFV case, this can also happen for lower values of 
82^^- In the same way, the enhancement of the total cross section through the new contributions from ct{tc) final 



states can be sufficient to achieve il^oh^ 



0.11. 



For completeness, we show in Fig. [7] the masses of the two lightest up-type squarks, the gluino, and the lightest 
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FIG. 7: Masses of the two lightest up-type squarks, gluino, and hghtest neutraUno (left) and flavour decomposition of lightest 
up-type squark (right)as a function of ^2^^^ for mo — 1500 GeV, mi^2 = 680 GeV, Ag = —500 GeV, tan/3 — 10, and fi > 0. 

neutralino as a function of the NMFV-parameter ^23^^ well as the flavour decomposition for the same scenario as 
discussed above. The squark mass splitting is increased due to the additional off-diagonal entries in the mass matrix, 
so that the mass of ui decreases. For large flavour mixing, it comes close to the neutralino mass, leading to the 
important coannihilation as seen in Fig.|4l The masses of U2 {— cl), the neutralino and the gluino remain practically 
unaffected by the considered generation mixing. 

B. Non-universal gaugino masses 

When considering 5*0(10) grand unification theories (GUT), the properties of the SUSY breaking mechanism are 
related to the breaking of an SU{5) subgroup into the standard model gauge group SU{3) x SU{2) x U{1). The 
relations between the gaugino masses Mi {i = 1, 2, 3) at the unification scale are given by the embedding coefficients 
of the standard model groups in SU{5). In particular, the unification constraint Mi = mi/2 of the CMSSM can be 
relaxed without spoiling the unification of the gauge couplings. Three independent parameters are then needed to fully 
parameterize the gaugino sector. A possible set is the wino mass M2 together with the two dimensionless variables 
xi — M1/M2 and X3 — M3/M2. The case xi — X3 — 1 corresponds to the CMSSM discussed above. Previous studies 
have shown that non- universal gaugino mass models have an interesting dark matter phenomenology [3^ [35| . 

We start by showing the relevant constraints in the (mo, M2) plane for Aq = GeV, tan/3 = 10, > 0, xi = 1/2 
and 2:3 = 7/4 in the upper panels of Fig. |51 The flavour violating parameter ^23 is set to zero (left) and to 0.95 
(right). In both cases, one WMAP-favoured region is situated around the resonance of the lightest Higgs-boson (for 
M2 « 300 GeV). The neutralino-stau coannihilation region is also present in both cases for mo w 100 GeV, next to 
the stau-LSP region. For M2 ^ 700 GeV, due to the large mass splitting for = 0.95, the squark-LSP region 

and its neighbouring coannhilation region are present, as it is the case in the CMSSM (see Fig. [3] right). In this case, 
a region with sizeable relative contribution from neutralino pair annihilation into tc is also present. Note that, due 
to the non-universality, the neutralino can be lighter as compared to the CMSSM. Therefore, this region is bounded 
at a certain value of M2, since below this bound the neutralino is not heavy enough to kinematically allow the top 
production. An upper bound for this region at a higher value of M2 is also observable, since above this value the 
neutralino is heavy enough to produce top quark pairs. As a result, this region lies in the range 450 ^ M2 ^ 850 
GeV. The region excluded by the 6 — >■ 57 branching ratio is sig nificantly larger for ^33^^ = 0.95 (mo < 1200 GeV, 
M2 < 150 GeV) than for 523^^ = (mo < 450 GeV, M2 < 100 GeV). For ^33^^ = 0.95 this region is similar to the 
one excluded by the light Higgs mass. 

We then show in the lower panels of Fig. |8]the constraints in the (xi ,X3 ) plane for a particular point of the parameter 
space for ^33 = (left) and $2^, = 0.95 (right). Note that the chosen point, for ^23 = 0-95 (i.e right panel) and 
xi — 1/2, X3 — 7/4, corresponds to a point where the relic density lies in the WMAP interval (see upper right panel 
of Fig. [HI). For S^f-^ = different allowed regions are visible. One is the very low xi (w 0.2) region in which the 
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FIG. 8: Top: Constraints in the (mo, M2) plane for xi = 1/2, X3 = 7/4, ^3^^^ = (left) and S^i = 0.95 (right) in NUGM. 
Bottom: Constraints in the {xi, X3) plane for mo = 320 GeV, M2 = 700 GeV, S^f^^ = (left) and S^f'^' = 0.95 (right). 



neutralino annihilates mainly to pairs of b quarks or tau leptons via a light Higgs resonance. Then there is a diagonal 
line corresponding to the heavy Higgs resonance (the neutralino mass increases with xi and the heavy Higgs mass 
with X3). For large xi (^ 1.8) the neutralino becomes mainly wino and annihilates strongly into W boson pairs, thus 
decreasing the relic density below the lower limit. In the favoured ellipse-shaped region for xi ~ 1.6, 0:3 ~ 0.8, the 
heavy Higgs resonance and W boson final states processes contribute in a way that the relic density is compatible with 
the required interval. For ^23^^ = 0-95 it is immediately visible that new allowed regions appear. In the low 0:3 region 
annihilation of neutralinos into tc (and tt for xi ^ 0.6) are the main contributions. However these new contributions 
essentially do not modify the shape of the allowed region compared to the ^23^^ ~ ^ case. Moreover, we can notice 
that this region is excluded by BR(6 57). The interesting region is for larger X3, where the neutralino coannihilation 
with the lightest squark appears. In this region the mass difference between the neutralino and the lightest squark 
is approximatively 30 GeV and the main contributions are neutralino-squark annihilation into top/charm quark and 
gluon, and squark annihilation into gluon pairs. While the neutralino mass increases with xi, the lightest squark 
mass increases with 0:3 until some point, and then decreases, which explains the particular arc shape of this region. 
Therefore this coannihilation region is relatively symmetric with respect to 0:3 « 1.1, and for each favoured point with 
X3 < 1.1 there is a corresponding favoured point with a;3 > 1.1. As explained below, the X3 < 1.1 part of this region 
is more strongly constrained by BR(6 — > 57), which explains our choice for the 2:3 value in the upper plots. 

Compared to the CMSSM, the predicted value for BR(6 — >■ 57) is much less constraining. First, because of the 
chosen value of Aq, but also because of the spectrum modification induced by the non- universal gaugino masses. This 
is illustrated in the lower right panel of Fig. [51 where mo and M2 have been fixed and the (dis)favoured regions are 
shown in the (xijX^) plane. It is striking that BR(6 -> 57) has a strong dependence on X3, which influences the 
whole SUSY spectrum through the renormalization group running. Most of the non standard model contribution to 
BR(6 57) comes from chargino and charged Higgs loops, the former being negative, and the latter positive. At 
low 2:3, the (negative) chargino contribution is dominant, which leads to a branching ratio far below the standard 
model prediction. With increasing 2:3 sparticles become heavier, and all absolute values of contributions decrease. 
However the absolute value of the chargino contribution come closer to the one from Higgs bosons, which leads to 
important cancellations and the branching ratio gets very close to the standard model value. For large X3, the chargino 
contribution (in absolute value) becomes smaller than the Higgs contribution. The sum is then positive, but rather 
small as the masses are large. 
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Among the other observables given in Tab. HI the ones which exclude some regions of the parameter space presented 
here are Aa^ and Ap, but only in the very low mass region which is already partially excluded by BR(6 57). 
For the plots in the (toq, M2) plane the neutralino mass excludes the region M2 < 250 GeV and the chargino mass 
excludes the region M2 < 150 GeV. The stop mass excludes also some region for which it is the LSP (i.e. already 
excluded). For the plots in the {xi, X3) plane, the neutralino mass excludes the region xi < 0.1. 

C. Non-universal Higgs masses 

Similarly to the mechanism leading to non-universal gaugino masses in 50(10) SUSY GUTs, depending on the 
exact representation to which the Higgs doublets belong, their corresponding SUSY breaking masses ttiho and ttihu 
need not necessarily be the same. In non-universal Higgs mass models they can therefore be treated as independent 
parameters at the high scale [H, H^l • For niHu = tuho = '^o the standard CMSSM is recovered. 

We start by studying the (toq, ^11/2) plane for fixed niHu — 1250 GeV, niHo — 2290 GeV, Aq = 0, tan/3 — 10, 
and /X > 0. The resulting excluded and cosmologically favoured regions are shown in Fig. [S] (upper panels), again for 
both the case of MFV (Jjs = 0) and a rather large flavour mixing parameter = 0.95. For = 0, the only 

allowed regions arise from the coannihilation of the neutralino with the superpartners of the tau or neutrinos and 
from the annihilation of neutralino pairs into VF^-bosons and top quark pairs due to the high higgsino component. 
The latter region is actually divided into two parts, parallel to the excluded region where certain squared sfermion 
masses become negative. 

For ^^3"" = 0.95, new contributions from XiXi ~^ ct{tc) make their appearance, as it is also the case in the discussed 
CMSSM and NUGM scenarios. Here, however, important coannihilations of the neutralino with the lightest squark 
are present, leading to a completely modified picture with respect to the MFV case. Among the WMAP favoured 
regions, only the one due to coannihilation survives the 6 — ^ 57 constraint. The discussion of this constraint (as all 
the constraints given in Tab. HI is here similar to the CMSSM one. 

Let us now study the {mnu, '^Hd) plane for fixed values of mo — 900 GeV and rai/2 = 700 GeV, shown in Fig. [S] 
(lower panels). For ^23^^ ~ ^ the situation is quite similar as described above. Two parallel allowed regions, where the 
neutralino is strongly higgsino and annihilating through the light Higgs resonance, are present. In addition, two other 
allowed regions, corresponding to the heavy neutral Higgs resonance, are observed for large M^u- Allowing for flavour 
violation {62^^ = 0.95), two large additional allowed regions appear, where neutralino-squark coannihilation processes 
are dominant (up to 50% including light squark annihilation into gluons). Note that the corresponding WMAP- 
favoured areas are very large as compared to =0. The flavour violating annihilation channel XiXi ~^ ciitc) is 

less important in this case. 

Again, all constraints have been checked to be fulfilled in the shown parameter space, except for the light neu- 
tralino/chargino mass in the region which lies close to the unphysical region. 

V. LHC PHENOMENOLOGY 

Finally, we discuss the collider phenomenology corresponding to CMSSM scenarios that feature new annihilation 
or coannihilation channels induced through flavour violating elements. Typical signatures for quark flavour violation 
in the context of squark production at hadron colliders have been discussed in Refs. [^-Q. A particularly promising 
process is the production of the lightest squark-antisquark pair, and their subsequent decay into charm- and top- 
quarks. The rather clean signature pp — > uiu\ — > cf(tc)xiXi, where the neutralinos would manifest as sizeable 
missing energy, might lead to up to 10** events at the LHC with ^/s = 14 TeV and an integrated luminosity of 100 

fb-i 

Alternatively, new contributions to the decay t ^ cj can increase its branching ratio to as much as 10^^ and thus 
render it detectable, e.g., in t< production at the LHC, as has been pointed out in Refs. [H, H^. In the case in which 
^23"^ is strongly constrained, e.g. by the Bg mixing, and only (^23^^ large, several years of high- luminosity operation 
might, however, be required. 

In order to evaluate the production of squarks and gluinos at the LHC, we have computed the relevant cross sections 
using the Monte-Carlo package WHIZARD 1.95 together with the associated matrix element generator O'MEGA [40l |. 
where the MSSM with the most general generation mixing as discussed in Sec. HI] has been implemented Q- We have 
employed the CTEQ6L [41] set for the parton distribution functions, the factorization scale being set to the average of 
the produced masses. Finally, the branching ratios of squarks and gluinos have been obtained using SPheno [l^ . 

In the first graph of Fig. [TOl we show the obtained dominant production modes of squarks and gluinos at the LHC 
with y/s = 14 TeV for the example scenario already discussed in Sec. IIV Al In the case of MFV {623^^ = 0), gluino 
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FIG. 9: Top: Constraints in the (mo, mi/2) plane for (Jj-j = (left) and = 0.95 (right) in NUHM. Bottom: Constraints 

in the [niHu, rriHo) plane for ^23™ = (left) and ^23^^ = 0.95 (right). 




FIG. 10: Dominant production cross-sections of up-type squarks and gluinos (left) and resulting NMFV signal cross-section 
(right) for the CMSSM scenario mo = 1500 GeV, mi/2 = 680 GeV, Ao = -500 GeV, tan /3 = 10, and fj, > as a function of 
the NMFV-parameter ^23'^''^ ^^e LHC with y/s = 14 TeV. The notation in the legend of the right panel is according to Eqs. 
(|26|l and (|28p . taking into account all possible combinations of (s)quarks and anti(s)quarks. 

pair production is dominant due to the colour structure, while squark-antisquark production is the subdominant 
channel. Due to the lighter mass, production of iii is preferred over the production of U2- We focus therefore on the 
production and decay of the lightest squark and the gluino. 

Due to only flavour-diagonal couplings, the gluino pair production cross section is practically independent of the 
NMFV-parameter 62^^- Contrary, the squark-antisquark pair production receives new contributions in a similar 
way as the neutralino pair annihilation discussed in Sec. IIIIl Since ui now has a sizeable c-admixture, initial states 
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containing charm-quarks can now contribute. Moreover, the hghter squark in the i-channel propagator enhances this 
channel. Finally, also the phase space is increased due to the decreased squark mass. Taking into account all these 
effects, the production of uiu* becomes the dominant channel for (^23^^ ^ 0.5 and reaches production cross sections 
of up to 10^ fb for S^f 0.95. For lower mixing parameters, gluino pair production remains numerically most 
important. 

Similar arguments hold for the associated production of a gluino and a squark. For the same reasons as given above, 
the production of uig is enhanced for large flavour mixing as compared to the MFV case. Note that, although the 
charge conjugated channel — ?> gul is not shown in Fig. [10] (left), it is taken into account in the following calculation 
of event rates. The associated production of g and heavier squarks is negligible in this context. The practically only 
flavour content of the second lightest squark 112 is t^, such that its production remains insensitive to the discussed 
flavour mixing in the right-right sector. 

The branching ratios of squarks and gluinos are also affected by flavour-violating elements, as can be seen from Fig. 
[TT] Since the lightest squark ui is a pure stop-like state in the case of MFV, it dominantly decays into ^ and 
{i — 1,2,3,4). Note that the decay into X2 is relatively small, since X2 is almost purely wino-like and couples only 
to the tiny iL-componcnt of ui. For increasing S^^ , the c^-content increases, and decays into final states including 
second generation quarks open up. Decays into heavier neutralinos are kinematically forbidden for 823^^ ^ 0.5. 
Although iji remains the dominant flavour in ui (see Fig. [7]), its decay into xft is closed for ^23^^ ^ O-^' since the 
mass difference between squark and neutralino is smaller than the top mass (see also Fig. [7]) . The only remaining 
decay mode is then ui — 

As a consequence, the NMFV signature 

pp^uiul^ct{tc)x°xi (26) 

discussed in Ref. cannot be realized for large 62^^ > 0.9. This channel can have sizeable event rates only for 
smaller flavour mixing parameters. This can be seen in the second graph of Fig. 1101 where we show the production 
cross section combined with the relevant branching ratios in order to estimate the signal rate at the LHC. 

For such important flavour mixing, a moderately sizeable signal rate can, however, stem from gluino pair and from 
gluino-squark production. In the former case, NMFV final states can be achieved through 

PP ^ 99 ^ cwi tui — > ccct XiXi, (27) 

where one of the produced gluinos decays into a top quark. As can be seen in Fig. Ill[ this decay mode remains allowed 
(and even dominant) for all values of ^23^^ . The notation of the final state in Eq. (P7|) is understood to include all 
possible combinations of quarks and antiquarks. 

The second possibility, mediated through associated production of a squark and a gluino, gives rise to signal events 
of the type 

PP uig cxi tui cct XiXi, (28) 
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where again the gluino decays into a top quark and the process is understood to include all possible combinations of 
(s)quarks and anti(s)quarks. 

In the second graph of Fig. [TOl we show the mentioned signal cross sections as a function of the NMFV-parameter 
'^23^^- discussed above, the signature in Eq. increases with 62^^, but drops when rriu^ — m^o < mtop, i-e. 
in the region where coannihilations with the lightest squark are most important. In this region, flavour violating 
signatures can be expected from the processes in Eqs. (I?f)) and (^5]) . They feature, however, cross sections that are 
smaller by about two orders of magnitude. 

For the LHC with y/s — 14 TeV and an integrated luminosity of 100 fb^^, the strongest signal in Eq. ((26|) can in our 
example scenario lead to up to about 10"* events. As discussed in Ref. 7], this signature is rather clean and not subject 
to important backgrounds. The region, where the correct relic density is achieved through efScient coannihilation, 
however, forbids this particular channel. The other potentially interesting channels may then lead to about a few 
hundred events each. Note that the event rate, the exact dependence on the NMFV-parameters, as well as the dark 
matter relic density are scenario-dependent. We can, however, expect a similar behaviour for other scenarios in the 
NMFV-MSSM. 



VI. CONCLUSIONS 

While the Minimal Supersymmetric Standard Model (MSSM) with a most general flavour structure has been 
extensively studied in the context of collider signatures, the possibility of squark flavour mixing has not been considered 
for observables related to dark matter so far. However, as the LHC is running and more precise cosmological and 
astrophysical experiments are taking data or being set up, it will become more and more important to take into 
account such effects when studying the interplay between collider and astroparticle phenomenology. 

In the case of neutralino dark matter in supersymmetric theories, flavour violating couplings can influence the 
(co) annihilation cross section, and in consequence the predicted relic density, in different ways. The strongest effect 
is due to the modifled mass spectrum of squarks, the lightest squark becoming lighter with increasing flavour non- 
diagonal terms in the mass matrices. The exchange of squarks in neutralino pair annihilation as well as the presence of 
coannihilation with a squarks become then important. Another effect comes from the fact that couplings of neutralinos 
to squarks are not diagonal in flavour space any more. This opens new (co) annihilation channels, such as XiXi ~^ 
or Xi"*^! ~^ ch^{cg,cZ^), which can give sizeable contributions to the annihilation cross-section already for moderate 
flavour violation parameters. 

Considering flavour mixing in the sector of right-handed up- type squarks, we have shown that the modified squark 
masses and flavour contents have a strong impact on the (co) annihilation modes. New annihilation channels are opened 
due to the presence of non-diagonal couplings in flavour space. These new contributions may become numerically 
important in particular regions of the parameter space. As a consequence, new regions that are compatible with the 
relic density constraint are opened. We emphasize the fact that these new regions are not excluded by the rather 
strong constraints imposed by flavour physics observables. Moreover, effects of lepton flavour violation on neutralino 
dark matter have recently been discussed in Ref. [Toj . 

A brief study of the corresponding LHC phenomenology has shown that the clean signature pp — ?> ctEl^^^^, that 
has recently been studied in Ref. |7|, can only be realized for moderate flavour mixing, when the lightest squark 
mass comes not too close to the neutralino mass. For rather large flavour violation, however, this channel is closed 
and NMFV-signatures arise through production and decay of gluinos rather than squarks. Such signatures include 
production of a top quark in association with charm-jets and may yield a few hundred events at the LHC with ^/s — 14 
TeV and an integrated luminosity of 100 fb^-'^. 

Since the annihilation cross section of the neutralino also governs the particle fluxes, flavour violating couplings 
would also have an impact on indirect detection of dark matter. In particular, additional c-~t mixing, as discussed in 
this paper, would change the spectrum of photons originating from dark matter annihilation. The impact of flavour 
mixing is, however, expected to be very small compared to the astrophysical uncertainties in this context. 

Direct dark matter detection might also be influenced by the discussed flavour mixing. Here, the scattering of a 
neutralino off a nucleus can proceed through squark-exchange, such that the charm-content in the nucleon becomes 
relevant if the lightest squark is a mixture of stop and scharm. In the same way, flavour mixing in the sector of 
down-type squarks would increase the importance of the strange quark in the nucleus. Detailed studies of direct or 
indirect detection of dark matter in the context of flavour violation are, however, beyond the scope of this work. 
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